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EM 

'  SS  =  RL(t)  +  TL(t)  -  ESL  -  101ogl0(7iCR(t)) 

•  SS  is  empirical  scattering  strength. 

•  RL  is  time- dependent  received  level. 

•  TL  is  time-dependent  two-way  transmission  loss. 

•  ESL  is  energy  source  level. 

•  7rCR(t)  ~  A(t)  -  insonified  area,  R(t)  is  range,  and  omnidirectional  source  is  assumed. 

'  Above,  terms  explicitly  dependent  on  time  are  highlighted  with  red  color. 

'  For  range-independent  environments,  SS=SS(0-j^^,0g^^j),  and  not  SS(t). 

'  Initial  approach: 

•  Use  TL  as  modeled  by  ARL/UT  based  on 

•  Inversion  of  environmental  properties  based  on  measured  38  g  shots  propagation. 

•  Forward  propagation  modeling 

•  Use  use  Wakeley  model  for  ESL  of  the  explosive  1  kg  source. 

•  Use  RL  as  recorded  by  VLA: 

•  Individual  sensor  for  angle- aggregate  SS. 

•  Beamform  to  extract  angular  dependence. 

'  Turns  out,  we  have  problems  with  source  model. 

'  Adjust  approach: 

•  Do  the  best  possible  under  circumstances  estimate  of  the  SS. 

•  Characterize  uncertainty  in  SS  resulting  from  uncertainty  in  source  and  TL. 


AV^rlz^  Qorr^n^ 


Scattering  Strength  Extraction 
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Source  Function  Sensitivity  to  Weight  Uncertainty 


Direct  VLA  data,  38  g  Shotl 

0r\ - 1 - 1 - A - ; - 


^V- 


- [■■^i  fi  I  -fc 


■I  ■.  V: - - 

i':.- 1.■:■■>^*^--■>^  ■;■■■: -r^r-  ■  ■  -■■ 


:.r  . ; :l} 


■-  ■ 'iV  ■ 


L,Vh - - 


. :-:■■  •  ; ' ' '  ■;•: 

. ’ .  . . 


3,6 


3.66 


3J 

t-TCn  S 


3J6 


3.0 


Source  PSD  depends  on  WBS  depth/mass. 

Depth  can  be  inferred  from  direct  arrival  structure. 

If  available,  pulse  separation  gives  estimate  of  weight. 

If  can  not  resolve,  assume  some  spread  of  weight, 
estimate  uncertainty  in  energy  density  spectrum  level 
(EDSL). 

Bottom  right:  EDSL  uncertainty  for  30%  weight 
spread. 


Estimi 

ite  depth 

Wakeley  model  -  estimate  weight 

Wakeley  model  -  quantify  uncertainty 


Uncertainty,  known  depth 


eiTor  range 
peridentO:! 


Eor  30%  weight 
uncertainty,  source 
level  uncertainty  is 
less  than  2.5  dB. 
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1  kg  Shot,  Direct  Arrival  on  VLA  and  Desensitized  Phone 


'  For  reference,  model  results  are  shown  for  1  kg  (red) 
and  1.8  kg  (green)  source  at  30  m  depth. 

'  Top:  desensitized  phone  data. 

'  Bottom:  direct  arrival  structure. 

•  About  30  m  detonation  depth. 

•  Individual  arrivals  are  hard  to  identify. 

'  Features  in  overloaded  direct  arrival  vaguely  hint  at 
1.8  kg/30  m  source.  But  1.8  kg  actual  weight  is 
unlikely  for  1  kg  nominal  charge. 

Can  infer  depth  but  not  weight.  I 

Use  nominal  weight,  quantify  error.  I 


1C 

20 

30 

4C 

m 

eo 

7C 

m 

90 


T - r 


[i  ■  v  —  .-Tu  Vi  .-^i  ■  ..T  ■'  ■  ■ .  -■  yi-.'iii  ^  .  f  ^ 

■X  -  -  !■  I.  W  H  yyur  -  U.-  I~1  .h-.l  ■  I  f  T.y,.-^ 


- IT  U  '.Vi-TV  W  U  -I  J’ 1  u '  jv'  .■.■  JUkA 

r . t  ^  .!-■  .V  --L--  -I*  4’  r-.  li," 

j  -iLI'  I  '.F^  linr.'  M.' ■■  '  I  ^  '  I  '  '  •  "i  ■  I  .  I  ■■I 

L. ...... .  .-r* ■*■.  rr.V"! vr'f v.v^.v^. a. -fii 

•.  ■.■■•>]■..■  .'J’ I  A'..:  t.  wjWKy  ,1  L.'  ..T-ll  VI.' I.  .  7  ■■  ■  _^-  ||■.■., 

r  i  ^  P 

;  J .:. J.  ' .Uh'  .  .■■I'lrTn.J  U^f  '.I'  A  {■.■.. 

^  M  hiiiV*>Art*rn  T/M^.  "Or  -i '  K  -"I  r  J 

'.7  rire  7ji  T:^  7.9  7J5  8 


vert_1l4gi_1.tsm 


BfJN  TECHNOLOGIES 


5 


AV^rlz^  Qon^n^ 


Range  Dependence  in  Source  Function 


1  kg  Charge 


22  Aiqai 


F,  Hz 
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Wakeley  model  to 
compute  source 
function  for  38  g  an(|i 
1  kg  source  at  50  m 
depth  at  range  of  lOp 
m,  1  km,  and  10  km 


Top:  EDSL,  2  Hz  arid 
20  Hz  resolution. 


Bottom:  EDSL 
deviation  from  that 
100  m. 

Left:  1  kg  charge. 

'  Right:  38  g  charge. 


200 


38  g  Charge 


10  10  10 
F,  Hz 

?7  Ain-m?  stfmfirnraRatqaDaMutenajTi  ESL  PSD  wawlev^a]  SQiiJiK 


at 


10  10 

F,  Hz 
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With  1  kg  charge,  observe  f-dependence  with  range. 

True  feature  of  nonlinear  source  or  Wakeley  model  artifact?  Implications  on  SS? 
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Single  Sensor  SS,  Mid- Water  Sensor  -  45  m 


Linear  TL  model: 

•  dependence. 

•  SS  increasing 
with  range  for 
420  Hz. 

Non-linear  TL: 

•  dependence. 

•  SS  decreases 
with  R  for  420. 

L/NL 

•  In  between:  F^-^. 

•  R-independent. 

Different  TL 

models: 

•  dB -difference 
not  huge,  but 
predict  different 
F-dependence. 

•  Which  F- 
dependence is 
right? 


BfJN  TECHNOLOGIES 


7 


AV^rlz^  Qon^n^ 


Single  Sensor  SS:  Range  Dependence,  Mid-Water  -  44  m 


Linear  TL: 

•  SS  in  dependent 
of  range  for  120, 
240  Hz. 

•  SS  increasing 
with  range  for  420, 
840  Hz. 

•  Increases  more 
for  840  Hz. 

L/NL  TL: 

•  SS  independent  of 
range. 

Is  it  an  accidental  fit, 
or  do  we  really  need 
to  consider  NL 
effects? 
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Question 


'  Assuming  non-linear  source,  can  “straighten”  range-dependence  of  SS. 

'  Does  this  mean  that  we  should  account  for  source  nonlinearity? 

'  May  be,  but  may  be  not,  as  there  are  potential  mechanisms  for  range  dependence  of  the  SS: 

•  Range-dependence  of  environment 

•  Scattering  into  different  angles  (coupling  between  modes)  increasingly  important  with  range 
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Single  Sensor  SS,  Shallow  Sensor  -15  m 


Results  very  similar 
to  mid-water  sensor 


BfJN  TECHNOLOGIES 


10 


AV^rlz^  Qon^n^ 


Single  Sensor  SS,  Deep  Sensor  -  70  m 


Results  very  similar 
to  mid- water  and 
shallow  sensors. 
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Uncertainty  in  SS  Estimate? 


'  Shot- to- shot  stability: 

'  Look  at  all  4  shots  provided. 

•  Shallow  sensor. 

•  Common  TL  model 

•  Linear  TL 

•  Range-independent  src. 

'  Very  similar  results  for  4  shots 


F,  Hz 


Observe  same  SS  for  shots  1-4. 

Potentially  2.5  dB  uncertainty  due  to  weight. 
Uncertainty  in  level  and  F-dependence. 

•  Linear  TL  model  results  in 

•  Higher  SS  above  200  Hz 

•  frequency  dependence  of  SS 

•  SS  increasing  with  R  for  R=4-10  km,  F>240  Hz. 

•  Nonlinear  TL  model  results  in 

•  Lower  SS  above  200  Hz 

•  F^-^  frequency  dependence  of  SS 

•  SS  independent  of  range  for  R=4- 10  km. 

What  is  the  correct  TL  model? 

Do  we  need  to  address  non-linear  souree  effeets? 
Conelude: 

•  dB -difference  between  L  and  NL  models  not  huge  - 
within  5dB. 

•  But  TL  model  affects  F-dependence  (important  for 
understanding  physics/mechanisms  of  bottom 
interaction) 

•  It  is  important  to  settle  on  the  correct  TL  model  and 
determine  correct  F-dependence  of  the  SS. 
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Array  Processing 


VLA  can  be  used  to  study  dependence  of  bottom  SS  on  the  scattering  angle. 

I  focus  on  bottom  subarray  (middle  column).  Reasons: 

•  Down  refracting  so  to  look  at  low  grazing  angles  need  to  be  close  to  bottom. 

•  Remove  contribution  from  surface,  focus  on  bottom  scattering  contribution. 

•  SNR  is  good  to  5-10  km  depending  on  steering  angle  and  frequency  band. 


Bottom 


250  Hz  octave 
Normalized 
by  horizontal. 

.Down 


"  j  Signal  clipped 
if  SNR<3  dB 


[,Pown 
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Angular  Dependence  of  Scattering  Strength:  Procedure 


'  To  determine  bottom  scattering  strength  SS(F,R,  a-,  a^),  need  to  normalize  beam  output  by 

•  Source  function 

•  Insonified  area 

•  TL 

'  Incidence  angle  a-: 

•  Can  not  be  controlled  with  explosive  source. 

•  Can  only  be  determined  from  modeling. 

•  May  change  as  function  of  range  due  to  stripping  of  high  angles. 

'  Scattering  angle  a^: 

•  Can  be  related  to  steering  angle. 


'  Using  different  array  beams,  can  obtain  SS(F,R  ,0C-Q,ag) 

•  otio  is  average  incidence  angle  (fixed  at  each  range). 

'  Caveat  1 :  As  for  single  sensor,  SS  potentially  depend  on  range  due  to  scattering  into 
different  angles  (coupling  between  modes). 

'  Caveat  2:  SS  may  depend  on  range  through  TL  dependent  on  angle  (bounce  loss).  Need  to 
use  TL(f,R,03),  where  0g  is  steering  angle.  Used  single-phone  TL(f,R)  instead. 


dB//m  I  P  dB/Zm 
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SS  in  Beam:  Range  Dependence 


240  Hz 


As  for  sensor,  have  similar  source-induced 
uncertainty. 

•  Linear  source:  dash  curves,  upward  rend 

•  Nonlinear  source:  dotted  curves,  downward  trend 

•  Linear/Nonlinear  combination:  solid  curves. 
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ss  vs.  Angle  at  Scatterer  and  Frequency 


SS(9o,f) 

•  f  is  frequency 

•  0G  is  average  array  beam  incidence  angle  at  the  bottom  (average  scattering  angle). 

Notation: 

•  Negative  angles:  results  observed  in  down-steered  beams. 

•  Positive  angles:  results  observed  in  up-steered  beams. 


-15  -ie 

Down 


^00 
Angle  Ht  Bottom:,  deg 


0  5 

Angli  at 


SS,  dB//m  I  p  SS,  dB//m 
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SS  in  Beam:  Angular  and  Frequency  Dependence 


Angular 

dependence: 

surprisingly 

Lambertian 


Very  steep  F- 
dependence: 

•  In  low  beams, 
steeper  than  F‘^ 

Error  in  SS 
extract? 

Or  tells  us  about 

scattering 

mechanisms? 
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SS  in  Beam:  Angular  and  Frequency  Dependence 


Lambertian  angular 
dependence. 

Steep  F- 
dependence: 

•  in  low  beams  - 
sub-bottom? 

•  Any  errors  in  SS 
extract? 

•  or  steeper? 
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Questions/Issues 


'  Assuming  non-linear  source,  can  “straighten”  range-dependence  of  SS  in  single  sensor  and 
in  beam. 

'  Does  this  mean  that  we  should  account  for  source  nonlinearity? 

'  May  be,  but  may  be  not,  as  there  are  potential  mechanisms  for  range  dependence  of  the  SS: 

•  Range-dependence  of  environment 

•  Scattering  into  different  angles  (coupling  between  modes)  increasingly  important  with  range 

•  Angular  dependence  of  TL  for  the  VLA  beam  data. 
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Conclusions 


'  Generally  low  bottom  SS,  increasing  with  frequency  for  f>400  Hz. 

'  Steep  frequency  dependence,  steepening  with  decreasing  of  grazing  scattering  angle. 

'  At  each  frequency,  angular  dependence  can  be  approximated  by  Lambert  Law.  Assuming 
15®  incidence  angle  (which  is  close  to  maximum  propagating  angle),  need  McKenzie 
coefficient  of  about  -50  dB  (at  120  Hz)  to  about  ^0  dB  (at  840  Hz). 

•  Low  values  of  McKenzie  coefficient. 

•  Lambert  Law  provides  functional  fit  but  does  not  correctly  describe  the  physics  of  scattering. 
Impossible  to  relate  required  values  of  McKenzie  coefficient  to  environmental  parameters. 


'  Need  to  compare  to  physics-based  theoretical  predictions  (e.g.:  surface,  sediment, 
volumetric,  SP/Born). 

•  F-dependence  of  SS  is  very  important  as  it  may  suggest  bottom  interaction  mechanisms 

•  For  example,  F^  can  be  argued  due  to  scattering  from  small  scatterers  distributed  near 
water/bottom  interface  (biologic  activity?) 

•  Current  uncertainty  in  TL  modeling  does  not  result  in  large  dB -difference  in  SS  extractions,  but 
affects  frequency  dependence. 

•  Need  to  better  understand  empirical  SS  for  better  understanding  of  scattering  mechanisms. 
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Importance  of  Nonlinear  Effects 


Pierce  ,  “Acoustics”:  Nonlinearity  counts  if  p-Po-27iFo/(p-C^)>a,  where: 

•  P=3.5  is  parameter  of  nonlinearity. 

•  Pq  and  Fq  are  pressure  and  frequency  characteristic  for  the  waveform. 

•  a  is  linear  attenuation  (in  nepers/m):  a^lO^/lOlog  10(e),  where  a  is  attenuation  in  dB/km. 

Use  empirical  expressions  from  Urick  for  a  charge  of  W  lbs  at  range  R  feet: 

•  Po=Co-2.16-104.(Wi^3/r  )i.i3  (Pa),  where  (0.454- 9.8)70.02542; 

•  Fq=10Vt  -  characteristic  frequency,  Hz,  where  t=58- (^wn/3/j-)0.22  _  characteristic  duration,  ps 
Use  empirical  curves  for  free  water  attenuation  from  Dyer: 

•  a=(3:30)  •  lO-^-  (Fq/1000)2  -  attenuation,  dB/km,  for  T=4^-22^  C  (lower  for  higher  temperature). 

•  a=a-  lO'2/lOlog  10(e). 

In  reality,  transition  to  linear  happens  earlier  due  to: 

•  a  and  a  are  higher  due  to  additional  attenuation  on  boundary  interaction. 

•  Pq  and  Fq  were  computed  with  account  to  spherical  spreading,  but  spreading  may  further  slow  steepening 
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38  g  Shot,  Direct  Arrival  on  VLA 


'  Small  source  ground  truth:  38  g  detonated  at  50  m. 

'  From  data:  detonation  depth  is  about  52  m  -  from  first  arrival  structure  on  VLA. 

'  On  direct  arrival  all  hydrophones  (including  desensitized)  clipped  but  multipath  and  bubble 
pulse  structure  is  clearly  evident. 

'  Separation  between  pulses  can  be  used  to  infer  detonation  strength. 


1-st  bubble 


Direct  WlA  data,  38  g  Shotl 
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38  g  Shot,  Direct  Arrival  on  Desensitized  Phone  vs.  Wakeley 


Below,  data  are  eompared  to  the  Wakeley  model. 

Levels  between  data  and  model  can’t  be  matched  because  of  clipping. 

To  fit  the  scale,  source  was  placed  at  100  m. 

•  Black  dot:  dB//|aPa  -  peak  level  at  100  m  per  Urick.  Shown  for  reference. 

•  Black:  data,  desensitized  hydrophone  (clipped). 

•  Red:  38  g  source  at  50  m  depth,  does  not  match  pulse  separation  in  data. 

•  Green:  25  g  source  at  50  m,  matches  pulse  separation. 

•  Blue:  38  g  source  at  60  m,  also  matches  pulse  separation,  but  depth  does  not  match  that  inverted 
from  first  arrival. 


Conclude: 

Actual  detonation  depth  is  about  nominal. 
Actual  charge  is  25  g  -  34%  less  than  nominal. 


